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Abstract 

The bound state problem for a gauge invariant quark-antiquark system is 
considered in the instantaneous rest frame. Focus here is on the long range 
non-perturbative interaction. A two-time Green's function is constructed for 
Salpeter amplitudes. The corresponding Schrodinger equation is found to be 
in the Salpeter form with a Wilson-loop term as the instantaneous kernel. 



1 Introduction 



The observed mass of a hadron results from the underlying dynamics of its con- 
stituent quarks and gluons. Understanding the process is made difficult however by 
the complexities of low-energy QCD which is widely believed to provide a mechanism 
for confinement. This is one of the most interesting and difficult studies in quantum 
field theory. The main problem stems from the essentially non-perturbative nature 
of confinement. This is so except in certain cases. For example, in the limit of heavy 
quark mass the theory is dynamically reduced to that of a single particle, and at 
short distances asymptotic freedom suggests that the binding is dominated by one 
gluon exchange which in the static limit yields a coulombic interaction. Thus heavy 
quarks forming a bound state remain largely insensitive to the details of confinement, 
and their static properties are well described by non-relativistic constituent-quark 
potential models Unfortunately, neither condition is met for the general case of 
deeply bound states, and a more fundamental approach is required. 
The starting point for a relativistic description of the meson spectrum is the covari- 
ant Bethe-Salpeter(BS) equation. It is the most orthodox framework for addressing 
the two-body problem in quantum field theory and has played a central role in the 
discussion for over three decades. On the other hand, the appearance of the rel- 
ative time variable lacks a clear interpretation in a Hamiltonian setting, and the 
usual approach is to make a three-dimensional reduction of the kernel while retain- 
ing relativistic kinematics. In addition, one or both single-particle propagators are 
often placed on an effective mass shell. Symmetric treatment of this last constraint 
leads to the Salpeter equation ^ whose basic statement on the analytic form of its 
single-particle components is given by ip^ = A^iIj± , also recognized as the no-pair 
condition 0. Apart from these issues are questions about the Lorentz structure of 
the confining kernel. Though the success of phenomenological models suggest and 
lattice simulations confirm static linear confinement, going beyond this limit requires 
knowledge of the kernel's dominant Lorentz components. But even this knowledge 
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would shed no light on specifics of its dynamical origin. 

What is needed and is essential for a discussion of this problem in non-perturbative 
dynamics is a manifestly gauge invariant formulation. This is provided by the path 
integral method. Defining the bound quark-antiquark amplitude in this way as 
a color singlet has led to some of the most fruitful work in the non-perturbative 
regime beginning with the pioneering observations of Wilson The simple picture 
offered in his area law for static quarks is compelling: Contributions to the meson 
propagator fall off exponentially in the area swept out by lines of chromoelectric flux 
joining the constituents world lines; hence widely separated paths are suppressed. 
The dominant effect is stated as an average over the transport of color along the 
contour of the minimum area, tlog{W{C)) ~ aSmin- The main difficulties here are 
technical and directly related to the fact that analytic expressions for the full single 
fermion propagator are not known. For example, completely static propagators lead 
to the well-known linear potential]^, and leading corrections yield 0(m^^) spin- 
dependent]^], and more recently as given in the works of Brambilla and Prosperi]^, 
spin- independent contributions. On another front, a full spinless propagator leads 
to relativistic flux-tube dynamics[§]. 

Here the following question is asked: Can the Lorentz and dynamical uncertainties 
present in the Salpeter approach be resolved and the technical difficulties of the path 
integral method overcome in a complimentary formulation which retains relativistic 
kinematics and spin structure? The answer seems to be in the affirmative with a 
resulting Hamiltonian similar to that of a relativistic fiux-tube model|P]. 

2 Two-body propagator 

The beginning parts of this section follow closely the treatment and notations found 
in ref[0 where the starting point is also the gauge invariant four-point Green's 
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function 

G{xi,X2,yi,y2) = ^(0|TV^2(a;2)f/(a;2,a;i)V'i(xi)V^i(?/i)t/(?/i,?/2)V^2(z/2)|0) (1) 
= ^{trU{x2,x,)Si{x,,yM)U{yi,y2)C-'S2{y2,X2\A)C). (2) 

The i/j 's here are fermion fields in the Heisenberg representation , U straight-hne 
path ordered exponentials, and 5* single particle propagators in the external gauge 
field A. C is the charge conjugation matrix while c denotes charge conjugation. The 
average is taken over gauge fields with virtual quark loops ignored. Positive and 
negative frequency components of the full single particle propagator 

iS{x,y\A) = e{xo~yo)S^{x,y\A) -e{yo- Xo)S'{x,y\A) (3) 

where S-^{x,y\A) = {0\^^+\x)^^~\y)\0) and S^{x,y\A) = (0|V5W(?/)^(-)(x)|0) , 
satisfy the homogeneous Dirac equation 

{iP,-m)S^{x,y\A) = (4) 

with Cauchy condition 

'S+ix,y\A) + S'ix,y\A)] = 7o5'(x-y). (5) 

L J xo=yo 

In the instantaneous Salpeter approximation ^^"^^ are expanded over free ± energy 
Dirac solutions, respectively, and thereby obey no-pair conditions which in the 
presence of the external field takes the form 

= A±(7r)^/'W (6) 

where tt = —I'V — gA, and A±(p) = {Eq±Hq)/{2Eo) are the usual energy projection 
operators with Eq = (p^ + m^)^/^ and Hq = 7*^(7 • p + m). The boundary condition 
(^) is then modified to 

S^{x,y\A)jX,=y, = A±(7r)53(x-y) (7) 
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so that each term in independently has the property of a propagator. By standard 
methods |TD| the solutions to satisfying (|^ are expressed as the phase-space path 
integrals 

S^{x,y\A) = rV[z,p]TAT,A±{p-gA) (8) 



X 



exp |«^^ dt[p ■ z T Eo{p - gA) - ^Ao]| 7° . 



In this expression Ta time orders gauge-field operators for xq > yo and antitime 
orders for xo < yo , while Tg does the same for Dirac matrices. Terms relevant to 
the minimal-area relation enter upon a semiclassical reduction. It is convenient to 
first isolate Wilson-loop factors by translating the integration variable p ^ p + gA. 
The momentum integration is then performed in the Gaussian approximation around 
stationary paths z± = ±p/i?o(p) yeilding 

Stiix,y\A) = J^V[z]TAT,A'^expi^tJ^^ dt[Tmil - z^y/^] - tg J^' A^dr^"^ ^'^ (9) 

The primed projector stands for evaluation at Pq = ±mz±(l — z^)"^/^ and an 
unimportant change in the integration measure has been suppressed [0. Consider 
now the action of the four-gradient 

id.S^i = ip,)dSti. (10) 

Equation (^) with constraint (j^) follow from the quantization of ([T0| ) by the obvious 
transformations, S^i ^ , pd — > — «V and Az^poS^ . Equation (|^) along 
with the relation C^^S^{y,x\A)C = [—S'^{x,y\ — A'^)Y , where r is the transpose 
operator, combine in the two-time Green's function to give 

G(xo;xi,x2|2/o;yi,y2) = 9ixo - yo) G+ + 0{yo - xq) G_ (11) 

with definition 

G±(xo;xi,x2|?/o;yi,y2) = (12) 



V[z,] P[z2]T,AfAf exp{z T dtY^[Tm,{l - z])''^] + \og{W {C^))}^^' 
yi •^y2 Jyo 
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for the closed contours C± 

yi 



xl 




+ 



X2 X2 

fig(l) 

time increasing to the right. Here the compatibihty of the instantaneous Salpeter 
picture with that of the Wilson loop is apparent: In addition to the intermediate 
propagation of two particles, both allow for the propagation of two antiparticles 
(double Z-graphs in time-ordered perturbation theory ) evidently represented here 
by . Single Z-graph components do not enter. As mentioned earlier, the straight- 



line approximation to the minimal-area law, tlog{W{C±)) = cxSn 
± dtLg is taken. In the center-of-momentum frame 



a\r\ / ds'j ^\ 
Jo 



(13) 



the exponential argument in ([l2|) can be written 



=F« / dt{Lq 
'yo 



=F« / dtL. 

'yo 



where 7(v) = (1 — v^) and = sxi_l + (1 — s)x2_l with f = (5ij — fifj)vj as the 
i-th component of v perpendicular to r = Xi — X2. With Lq = I]^=i rrijil — z|)^^^ 

(14) 

(15) 

(16) 
(17) 



From the stationary condition, action of the four-gradient yields 



. dxn 



lA- 



= J2 ±"^jXj±7(xj) + o-|r| 

j 



dsj(yt)^t 
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and 



r)r 

^^^,{^A^) = ±{T.^j±-^-L) (18) 



= ±{^mj-f{icj) + a\r\ f'ds^ivt)) (19) 

. JO 

= ±(Em,7(x,) + y). (20) 



Then 



j:-idjG^ ^ ($:P,, + P,)G± (21) 

j j 

id,,G± ^ {H;,, + H',^±V±)G± (22) 



3 Eigenvalue equation 

Given ri, r2) and ri, r2) as positive and negative frequency compo- 

nents of the single-time bound state qq amphtude with t < t', the freely propagating 
amplitude at intermediate times is found with help from the Green's function 

^r2[G{xo; xi, xsjt; ri, rs)^^^^^; ^i, r2) (23) 
-G{xo; xi, X2|i'; ri, r2)x^"-*(^'; ^i, r2)] 

and so satisfies the time-independent Schrodinger equation 

id,,X - Hx (24) 
where the Hamiltonian is given by 

H = H',, + H;,^ + V+-V.. (25) 
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Here Hq represents the quark's kinetic energy and V the gluon field contribution 
from the Wilson- loop (which apparently enters as the instantaneous BS kernel). 
Projection operators on V prevent mixing with the non-normalizable continuum 
states of the same energy. 

This result should be compared with others from the Wilson-loop formalism. The 
classical limit to the spinless Hamiltonian of ref|Q is recovered from (^) by the 
obvious reduction, Hq —>■ +Eq, which incorporates the no-backtracking constraint; 
that is, the second loop of fig(l) does not contribute ( — > 0) in this approxi- 
mation. 0(l/m^) terms of {H)i2 in the standard Dirac representation give leading 
relativistic corrections to the static ar interaction; the reduction is equivalent to a 
Foldy-Wouthuysen transformation. In the center of momentum frame the correc- 
tions are 

Korr ^ -a{l/ml + l/ml-l/{m^m2))Ly{6r) (26) 

-a{L ■ si/ml + L ■ S2/ml)/{Qr) . 

The spin-independent term above agrees with the semirelativistic result of ref[|^ 
omitting a angular-momentum independent contribution resulting from their par- 
ticular operator ordering prescription, Vc/(r, p) — > {Vci}ord = Vqm- An alternative 
prescription better suited to quantization of (pSj) would be in terms of classical co- 
ordinate and velocity observables, (r, v). In this case the v_|_ operators are found 
as truncated matrices in a suitable basis from the symmetrized orbital angular mo- 
mentum equation, L(r,v^) = x, x pj, for a given state. The energy eigenvalue 
equation thus becomes a matrix equation in the radial coordinate only and is solved 



variationally. For details of the method see e.g. ref|ll 



Spin corrections responsible for fine and hyperfine structure of the spectrum are of 
course an important and delicate concern. They must be handled carefully. The 
numerical coefficient of the spin-orbit term in (p6D , -1/6, is at variance with the 
-1/2 factor found in most of the literature. Most often the -1/2 follows from an as- 



sumption of dominant scalar confinement [12] which in relativized models based on 
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a BS reduction is thought to be spectroscopically favored over other Lorentz forms. 
On the other hand, this same disagreement appears with the 0{l/m?) Wilson- loop 
result of ref[0 where no such assumption is made. From a purely mathematical 
point of view this results from the different ways in which the holonomic straight- 
line condition is applied to the minimal- area relation; a factor of -1/6 (-1/2) results 
when the constraint is imposed before (after) the variation implied here in equa- 
tion(|?I|) and in appendix B of ref0. Only the former option, followed here, is in 
line with proper procedures for the mechanics of constrained systems. It should be 



pointed out that Gromes has deduced a relation [0 in support of the -1/2 factor 
(and scalar confinement) from arguments of Lorentz covariance following a reasoning 
in Representation Theory on the group structure of Poincare transformations. The 
relevance of this formalism in the present rest-frame Hamiltonian context is however 
not entirely clear and the question will not be entered into here. 



4 Summary 

Salpeter amplitudes have been introduced into the path integral formulation of QCD 
beginning from the gauge-invariant four-point function. A Hamiltonian description 
for mesons as bound quark-antiquark states has been derived in a systematic and 
straightforward manner. Interestingly, a time derivative of the Wilson-loop operator 
appears as the instantaneous BS kernel. Both relativistic kinematics and spinor 
structure of the amplitude have been preserved. The focus here has been on long 
range non-perturbative effects of the gluon dynamics, though a realistic calculation 
of the spectrum must take medium and short range contributions into account as 
well. These enter asymptotically as a coulombic interaction which is easily added 
to the present result. An apparent discrepancy between the 0{l/m?) spin-orbit 
coefficient and the relation of Gromes has been noted and is under study. 
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